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ABSTRACT 


A theoretical model is developed to establish an adhesion criterion 
for the interaction of a gas atom or molecule with a nonmetallic sur- 
face. As a first approximation a linear lattice mass points of the 
solid surface is considered and the classical lattice dynamical model 
is used. 
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INTRODUCTION 


12 3 

Classical as well as quantum mechanical approaches * ’ have been 
developed by different investigators in arriving at satisfactory 
models for the trapping of gas atoms by solid surfaces. In this in- 
vestigation the classical approach using the formalism of lattice 
dynamics is used. As a plausible potential for the interaction be- 
tween the surface and the impinging atom, the well-known Lennard- 
4 

Jones potential is used. 

The intermolecular interactions are represented by appropriate 
spring constants; a uniform spring constant for the binding between 
the atoms of the solid and a different value for the interaction be- 
tween the incoming gas atoms and the nearest surface atom are used 
in this approach. 

The required adhesion criterion is developed in the form of an 
estimated potential well depth which is used to compare the kinetic ' 
energy of the incident molecules. It is assumed that molecules with 
energies greater than the well depth will be reflected back whereas 
those with energies less than this amount will be absorbed. 

The one- dimensional semi-infinite lattice model used in the 
present discussion can be extended to the more realistic three- 


dimensional situation. 
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THEORETICAL APPROACH 
One dimensional chain lattice 

The solid surface can be approximated to a semi-infinite one 
dimensional lattice with identical point masses at each lattice 
site holding each other through a potential which can be represented 
by a suitable spring constant k. The interaction is limited to only 
the nearest neighbors. Any change in energy at any of the lattice 
points will be transmitted to nearest neighbor and from that to the 
next and so on. 
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Fig. 1. one dimensional lattice 


Let the instantaneous position of the n and (n + 1) atom 

in the chain be denoted by X and X respectively. Further, let 

n n + 1 

the equilibrium distance be denoted by d. 

Let there be a disturbance in the linear chain by an approaching 

impurity atom (gas atom) leading to a displacement <f> for the n c 

n 

th 

atom, thereby bringing the n atom to its instantaneous position de- 
noted by X . Under these circumstances one can write 
n 

X n = nd + <j> n . 1 

In a similar manner one can write for the (n + 1)^ atom 

X =(n + l)d + <f> . 

n + 1 n 


2 
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The interaction energy between the n t ^ 1 and (n + atom will 
be a function of the positions and can be represented by a potential 
function of the form 

V(x) - V(jX n + x - X n i). ' 3 


Equation 3 can be rewritten using equations 1 and 2 and one 
obtains 

V(x) = V |[(n + l)d + <J» n + x ] - [nd + <j» n ]| 

4 

■ V {d + !♦„+ 1 ' *nl> • 


The change in displacements between successive lattice mass 
points will be much smaller than the mean distance between successive 
lattice sites. 


U , , - <f> I « d 
,Y n + 1 v n> 
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Therefore, it is reasonable to make a Taylor expansion of the potential 
function v(x) about x=d which leads to 


V(x) = V(d) + |x=d (U n + 1 - <f>J) 

+ -^| • |(U n ' + q - 4> n l) 2 + 

3x x=d 

In equation 6, v(d) is a constant Since d is the equilibrium dis- 
tance, the potential function should be an extremum at x=d. This would 
mean that. 
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3v 

9x 


= 0. 

x=d 
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At equilibrium position force is zero. Using the analogy from 

spring, one can assume a relation of the form fa - kx or — akx. 

3x 

3 

This would imply that at — — i will be a constant. This simplifi- 

3x 2 x=d 

cation yields from equation 6. 


r w - c i + 4 k( K + 1 - *„i )2 • 


At the site of the n t41 atom, there is a similar contribution 
til 

from the (n-1) atom. Thus the total potential at the site of the 
n t41 atom can be written as 


>w - + j - +„l > 2 + - * n - il>" + c - 9 


The force on the n L lattice atom is then written as 


= - 21 


n 


rn 


- - K «n + 1 - ' +„ - 1> 
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Or 


■ K <*„ + 1 + *„ - 1 - 2 V 


The corresponding equation of motion is written as 
d 2 <j> 

m -—2 = K (<j> + * - 2* ) 

dt 2 n + 1 n-1 n 


11 
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Equation 11 is a representative equation which can be used to 
study the propogation of energy transfer through a lattice medium. 
The above considerations can be further simplified by adopting the 
following notations: 


Let to 2 -£•; 
0 ra 


12 


And T = w^t . 


13 


Then, -77 = ^ 


d^ _ d^ . dt _ 


dt dt dt 


0 3x * 


Similarly, 


d 2 <i> 

dt 2 


d ,d<h 
“0 dt W 


“o ft 3x <»> 


= 0 ) 


2 

0 dt 2 
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Using equation 14 in equation 11, one obtains for the (n + 1)^ 
particle 


mio r 


+ 1 ~ K(<f, n + 2 " 2 * n + 1 + V » 


and 


♦„ + 1 * *n + 2 ' 2 K + 1 + *n 


15 


Equation 15 represents the motion of the (h + l) t ^ 1 atom in terms of 
the displacement of the nearest neighbors. 
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When an external impurity atom approaches the zeroth atom, 

the motion of the coupled system can be written as given in 10 if 

the relative displacement of the incoming atom is denoted by and 

the spring constant representing the coupling between the zeroth 

atom and the trapped atom is denoted by K , then equation 10 leads 

o 

to , 

2 

m— 0 = K(<j) 1 - <fr Q ) + K 0 (£ - <{> 0 ). 16 


This is the equation of motion of the zeroth atom on the surface. 


d^<j) 

moiQ — 2“ = Kfo], " $ 0 ) + K q (£ - <J> 0 ) , 

dt 


K„ 


Or 


7 w o = ( *i ~ <i> o ) + r (? ~ V 


17 


K 2 

Since — is defined as u> , and if one further defines 
0) o’ 


K r 


K 


as g. 


equation 17 reduces to 


4* = ^l ” *o> + P _ 4> 0 ) 


18 


The motion of the impurity atom alone can be described by 

m 0 I = - K q ^0 “ ^ ’ 

dt . u 

— £ = K 0 m(<j> - K) , 

K 0 K~ 0 

m Q ? = Bm(<j> 0 - 0 . 19 


Or 

Or 
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If one writes m o /m = V, equation 19 reduces to 
y? = 6(4> 0 - O = By(r). 


One then rewrites equation 18 in terms of equation 20 to obtain 

♦ o = ^1 " + VK, 21 

y(r) = 4 >q(t) - 5(t) 

and y = ^' 0 - 22 

Equation 21 can be rewritten in terms of 20 as 

♦o m '^ 1 - V + 23 

Equation 23 when applied in equation 22 one obtains 

• y - (+! - V + PyCO - 5 (t) 

= (4*^ - <f>o + 3y(r) - y(r), 24 

y + 1 

Or y(r) = (4*2^ - 4>o^ “ &( — “ — ) y(f); 25 

r* 

Or y(r) + g( y(r) = (^ - 4 >q * 26 

Equation 26 represents coupled equation of motion. This is similar 
to a harmonic equation with a forcing term. At this stage it is 
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convenient to revert the variables back in terms of real time. 

—2 7I +8( T Ll >^'+l-*0- 

“0 dt 

2 2 M + 1 

Let « = Uq 3(— —-) 


and (4> 1 - <p 0 ) = f (t) 


Equation 27 is now rewritten as 


2 

d y , n 2 

+ ft y = 

dt / 


The factor which modifies the free lattice frequency is {g[(u+D/y] } 


and this is the effect of the coupled external molecule. Equation 28 is 
a harmonic equation with a forcing term and the forcing function is 
f (t) = Wq ( 4> 1 - (f> 0 ) • 


TRAPPING ENERGY 


Molecular potentials 


In a simple classical model, the molecular interaction potential 
can be expressed by various different potentials such as Lennard-Jones^ 
which is usually written as 

V (r) = 4e[(-f-) 12 - (~^) 6 ], 30 


where r is the distance from the surface lattice atom, a is the distance 


of closest approach by the incoming atom and e is the depth of the poten- 


tial. 
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A more realistic potential would be the Buckingham 3 potential 
which includes the induced dipole-induced dipole interaction and 
the induced dipole-induced quadrupole interactions. Buckingham po- 
tential is usually written as 


. , , . -ar -6 ' -8 

$(r) = be - cr - e r 


31 


While it is a more realistic potential than the Lennard-Jones (6-12) 
potential, it is more difficult to handle. A third potential which 
is flexible and often used in calculations is the Modified Buckingham 

g 

(6-Exp) potential which has the form 

Mr) = [-|exp( [1 - rsH) " 

1 a r> r 

— max 

$(r)=a- r << W 32 


This is somewhat more flexible than the Lennard-Jones potential since 
it permits the variation of the low velocity collision diameter. Another 
potential often encountered in literature is the Stockmayer Potential^ 
which is a superposition of the Lennard-Jones potential. 

q 12 o 6 ^a ^b 

Mr, e a , e b , <t> h - <j> a ) = 4s[(— ) - (— ) ] — g(e a ,e b ,<j> b ,-<f> a >- 33 

r 


The last term in the expression on the right-hand side of equation 33 
accounts for the angular dependence of the dipole-dipole interaction. 



Kr) 


- 10 - 


This potential describes well the interaction between such polar 
molecules for which dipole-quadrupole and higher multipole inter- 
actions are not important. 



Fig. 2. Spherically symmetric potentials 


The interaction of an adsorbed impurity atom by a solid surface, 
it is assumed, can be best described by the Lennard-Jones potential 
which is shown in Fig. 2b. The potential itself is described by 
equation 29. It can be reasonably assumed from the shape of the 
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potential, that the intercept 0 can be taken as the distance of closest 
approach for an external molecule. One further makes the stipulation 
that the depth of the potential well, e, can be interpreted as the 
mean binding energy of an incoming atom when it is absorbed. Under 
these circumstances, it is reasonable to assume that an incoming mole- 
cule will be trapped if it comes closer than the mean distance of 
closest approach and if the magnitude of its kinetic energy is less 
than the potential depth. This is a classical concept. It must be 
noted that the interaction potential does not become a trapping poten- 
tial until the approaching molecule is closer than Rg where Rg is 
shown in Fig. 2b. 

An estimate of trapping energy 

Equations 28 and 29 may now be used to estimate the trapping 
energy Eg as follows. A simplifying assumption that one can make is 
that the lattice response forcing function, which is given by 


,, \ 2 
f(t) = (Oq 


(4>i - 


is small compared to the relative acceleration 



and the average dis- 


o 

placement y. This would lead to further simplifying situation that 
the system would oscillate essentially with free angular frequency 0. 
Further, if one approximates equation 29 to be of harmonic form, which 
can be written in the form 


V' = 1/2ct 2 


34 
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Obviously one can then write the expression for the lower potential 
(Fig. 3) in terms of the harmonic potential as 

V = V’ - E . 36 

B 


Using equation 34 in 36 and evaluating the potential for the lower 
curve (Fig. 3) for r is equal to a, one obtains 


0 




37 


E 

B 




38 


One makes use of the relations given in equations 28 and 35 to reduce 
equation 38. Thus, 


2 

k = m^ 


2 u + 1 
m Q to [B(- )]. 


39 


From Fig. 3, 


n 


a 
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The use of equations 39 and 40 in equation 38 leads to. 


E B = \ m 0 -0 te(L T^ )](R B " a) ‘ 


41 


Obviously, equation 41 is an expression for the binding energy in terms 


of the quantities on the right hand side of equation; where m^ is the 

K 1/2 

mass of the trapped molecule, w n , which is equal to ( — ) , is the frequency 

v m 



term, K being the lattice spring constant, m the mass of lattice atom; 

3 is the ratio Kq/K where Kq is the equivalent spring constant between 

m 0 

the trapped atom and the surface atom; p or — is the ratio of the 

m 

masses of trapped and surface atoms; R the mean molecular trapping dis 

B 

tance and a is the distance of closest approach. 

It is obvious that if the quantities on the right hand side of 
equation 41 can be evaluated that would provide the necessary informa- 
tion on the probability of adsorption of an incoming molecule by a 
known target surface. 

Distance of closest approach 

In the preceding discussion several simplifying assumptions have 
been made in order to arrive at an expression for the trapping or bind- 
ing energy E as given in equation 41. The more important parameters 
B 

which are to be determined to evaluate E are o, the distance of 

B 

closest approach; R^, the mean molecular trapping distance; and Kq, 
the intermolecular spring constant between the trapped atom and the sur 
face atom. One makes further simplifying assumptions to evaluate these 
quantities. The quantity a which is the collision diameter can be 
reasonably defined to be the distance of closest approach between two 
colliding molecules. The value of r, occurring in the Lennard- 
Jones potential expression, for which <J>(r) becomes zero, will be a 
measure of o. The collision diameter is related to the molecular co- 
volume 'bQ 1 through the relation, 
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An alternate method to evaluate a is to use the diffusion data for 
molecules. However, the results obtained thus are in reasonable agree- 
ment with values obtained from the molecular co-volume as given in 
equation 42. 

In published literature there are available values of molecular 
co-volume and o of various molecules. Table 1 is adopted from Curtiss 
and Herschf elder to show some typical values of these quantities. 


Compound 

0 

(A° ) 

b 0 

cc/mole 

chci 3 

2.98 

33.45 

chci 2 f 

4.82 

141.0 

c 2 h 5 ci 

5.41 

199.7 

h 2 o 

2.65 

23.45 

nh 3 

2.60 

22.12 





Table 1 


Sample Values of a and b 


0 


In terms of the constants appearing in equation 42, one can express 
b Q = 1.26515a 3 . 


molecular co-volume as 
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The trapping distance is of interest here now. From Fig. 3, it 
is seen that when r is equal to Rg, 

V(r) = V(Rg) = -e = -E . 43 


This information can be used in the Lennard-Jones potential 


a 12 a 6 

V(r) = 4e[ (— ) - (— ) ] 

r r 


From equation 43, this becomes 


a 12 a 6 

-e = 4e [ ( ) - ( ) ]. 

R R 

B B 


44 


Equation 44 leads to quadratic of the form 


a 12 a 6 

( ) - (— ) + -f- - 0. 

R R 4 

B B 


45 


Or 


[( "T )6 " “f ]2 = °* 

r b 2 


This leads to two equal roots 


(-) 6 - A 

V R 2 

B 


46 
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This leads to a solution 


1/6 

\ = (2) o, 

or 

R b = 1.122a. 


47 


When the value of Rg from 47 is used in equation 41, one obtains 
for binding energy 


1 2 u+i 2 

E = — m w [0(— )] (.1220) . 

B 2 0 0 y 


48 


It remains to calculate the spring constants K between the lat- 
tice mass points of the target and Kq, the similar parameter relating 
the incident molecule to the surface atom. In the latter case, where 
Lennard-Jones potential is adopted, Curtiss and Herschf elder have 
elaborately discussed the applicability of thermal diffusion data in 
conjunction with ordinary diffusion data to determine a and e which 
occur in Lennard-Jones potential which describes the interaction of a 
pair of dissimilar molecules. On the other hand, this can be treated 
as a dipole-dipole interaction in which case the force of attraction 
can be written as 

y i y 2 


49 
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where and p 2 are the dipole moments and r the mean distance. Since 
we have replaced this force by a spring, equation 49 can be written as 


W = Vs 


50 


where R is the separation distance at which the attractive potential of 
s 

the surface atom is experienced by the incoming molecule. 


Equations 49 and 50 immediately lead to the result 


K o" 


6h]h2 


51 


Re 


PjL and induced dipole moments, can be evaluated from molecular 

data. The separation distance R can be approximated from the shape 

s 

of the potential in terms of q and a. From fig. 3 one can write 

0 


R 0 = 2q + a = 2R„ - 0 . 
s o £ 


Finally, the lattice spring constant can be obtained from the properties 
of the material. Thus for example the speed of sound in a material is 
given by 


v 



where y is Young's modulus, and p the density. Using the relation between 
to and v the frequency, one can then write 

v = Am, 
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where X is the wave length. In the linear approximation of a spring 
type force between the atoms, then one can write 

m = 



or 


K 


2 

mm 


- f v \ 2 

- »<— ) • 


Thus the various quantities appearing in equation 48 are evaluated to 
obtain an estimate of the trapping energy. 
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